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Abstract 

A  coefficient  inverse  problem  for  the  non-stationary  single-speed  transport  equation  for 
t  e  (0,  T)  with  the  lateral  boundary  data  and  initial  condition  at  t  =  0  is  considered.  Global 
uniqueness  result  is  obtained  via  the  method  of  Carleman  estimates. 


Introduction 

The  transport  equation  is  used  to  model  a  variety  of  processes  of  particle  transport,  such  as 
neutron  diffusion,  scattering  of  light  in  the  atmosphere,  propagation  of  y  -rays  in  scattering  media, 
etc.  (see,  e.g.,  the  book  of  Case  and  Zweifel  [6]).  Coefficient  inverse  problems  (CIPs)  for  the 
transport  equation  are  the  problems  of  determining  of  the  absorption  coefficient,  angular  density  of 
sources  or  scattering  indicatrix.  from  an  extra  boundary  data.  They  have  a  variety  of  applications  in 
optical  tomography,  theory  of  nuclear  reactors,  etc.  (see,  e.g.,  the  book  of  Anikonov,  Kovtanyuk  and 
Prokhorov  [  1  ]  and  references  therein).  This  paper  addresses  the  question  of  the  global  uniqueness  for 
a  CIP  for  the  non-stationary  single-speed  transport  equation  with  the  extra  lateral  boundary  data. 
Stability,  uniqueness  and  existence  results  and  references  to  such  results  for  CIPs  for  the  stationary 
transport  equation  can  be  found,  e.  e.g.,  in  [1]  and  in  the  book  of  Romanov  [25].  Uniqueness  and 
existence  results  for  CIPs  for  the  non-stationary  transport  equation  were  obtained  in  the  works  of 
Prilepko  and  Ivankov  [22],  [23]  and  [24].  The  results  of  [22]  and  [23]  were  obtained  for  special 
forms  of  the  unknown  coefficient  using  the  overdetermination  at  a  point.  Also,  uniqueness  and 
existence  results  were  obtained  for  an  inverse  problem,  where  complete  lateral  boundary  data  is  not 
present  but  both  initial  and  final  conditions  (at  t  =  T)  are  given;  see  [24] .  For  some  recent 
publications  on  overdetermined  inverse  problems  for  the  transport  equation  see  Tamasan  [27]  and 
Stefanov  [26].  A  derivation  of  the  transport  equation  for  the  non-stationary  case  can  be  found,  for 
example,  in  [6] . 

The  proof  of  the  main  result  of  this  paper  is  based  on  a  Carleman  estimate,  obtained  by  Klibanov 
and  Pamyatnykh  [17].  Traditionally,  Carleman  estimates  have  been  used  for  proofs  of  stability  and 
uniqueness  results  for  non-standard  Cauchy  problems  for  PDEs.  They  were  first  introduced  by 
Carleman  in  1939  [5],  also  see,  e.g.,  books  of  Hormander  [7],  Klibanov  and  Timonov  [T9]  and 
Lavrent’ev,  Romanov  and  Shishatskii  [21].  Bukhgeim  and  Klibanov  [4],  [13]  have  introduced  the 
tool  of  Carleman  estimates  in  the  field  of  CIPs  for  proofs  of  global  uniqueness  and  stability  results 
for  CIPs,  also,  see  Klibanov  [14],  [15],  and  Klibanov  and  Timonov  [19],  [20].  This  method  works  for 
CIPs  with  single  measurement  data  for  the  time  dependent  Partial  Differential  Equations  (PDEs),  as 
long  as  the  initial  condition  is  not  vanishing  and  the  Carleman  estimate  holds  for  the  corresponding 
differential  operator.  Recently,  Klibanov  and  Timonov  have  extended  the  original  idea  of  [4]  and 
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[13]  -  [15]  for  constructing  of  a  globally  convergent  numerical  method  for  CIPs,  including  the  case 
when  the  initial  condition  is  the  5-function;  see  [19]  for  details  and  more  references.  A  variety  of 
works  with  uniqueness  and  stability  results  for  coefficient  inverse  problems  were  published,  which 
are  based  on  the  idea  of  [4]  and  [13]-[  15],  see,  for  example,  Bellassoued  [2],  [3],  Imanuvilov  and 
Yamamoto  [9],  [10]  and  [1 1  ]  for  the  case  of  hyperbolic  equations.  The  Lipschitz  stability  result  for 
the  CIP  for  the  non-stationary  single-speed  transport  equation  for  t  e  (-T,  T)  with  the  complete 
lateral  boundary  data  and  initial  condition  at  t  =  0  followed  by  the  corresponding  uniqueness  result 
was  recently  obtained  by  the  authors  [18]. 

A  natural  question  arises:  Is  it  possible  to  consider  the  CIP  for  the  non-stationary  single-speed 
transport  equation  t  e  (0,  T)  rather  than  for  t  e  (-T,  T)  and  obtain  at  least  a  global  uniqueness 
result  for  this  case,  which  would  be  similar  to  [18]  ? 

Since  there  are  no  other  techniques  so  far,  except  of  one  of  [4],  [13]-[15],  [19],  which  would 
enable  one  to  prove  global  uniqueness  for  the  multi-dimensional  CIPs  with  the  single  measurement 
data,  we  have  no  choice  but  to  use  that  method.  However,  this  method  doesn't  necessarily  allow  one 
to  consider  the  problem  in  (0,  T)  with  the  initial  data  at  {t  =  0}.  The  obstacle  is  due  to  the  presence 
of  an  integral  in  the  Carleman  estimate,  depending  on  the  divergence  term  V,  where  the  integration  is 
carried  out  over  {t  =  0}  (see  Lemma  1).  In  the  case  of  hyperbolic  inverse  problems  this  integral  is 
zero  due  to  the  Carleman  estimate  for  the  principal  part  of  the  hyperbolic  operator,  see  Theorem 
2.2.4  in  [19].  This  is  basically  due  to  the  fact  that  the  hyperbolic  operator  contains  the  derivative  of 
the  even  order  (2)  with  respect  to  t.  Thus,  it  is  possible  to  consider  hyperbolic  inverse  problems  in 
(0,  T).  However,  in  the  Carleman  estimates  for  parabolic  and  transport  equations  such  integrals  are 
not  zero.  Thus,  these  Carleman  estimates  do  not  allow  to  consider  the  corresponding  inverse 
problems  in  (0, 1),  at  least  directly. 

Another  method  to  try  is  to  consider  the  forward  problem  in  (0,  T),  construct  an  extension  of  its 
solution  into  (-T,  0)  and  consider  the  inverse  problem  simultaneously  in  (0,  T)  and  (-7’,  0).  So  that 
the  integrals  over  {t  =  0},  arising  from  the  divergent  term  V  in  both  cases,  would  cancel  out.  In 
order  to  use  the  same  Carleman  estimate  for  the  problem  in  (0,  T )  and  (-7’,  0),  the  principal  part  of 
the  differential  operator  should  not  change  in  such  an  extension. 

It's  quite  easy  to  find  such  an  extension  for  the  hyperbolic  inverse  problem.  For  example, 
considering  the  equation 


utt  =  An  +  a(x)u,  n(x,0)  =f(x),  ut(x,  0)  =  0, 

we  see,  that  the  extension,  even  with  respect  to  t,  leaves  the  form  of  the  operator  unchanged  in 
{-T,  0).  However,  this  is  not  the  case  for  parabolic  problems.  Indeed,  considering  the  simple 
parabolic  equation 


ut  =  Au  +  a{x)u,  u(x,  0)  =  fix), 

and  trying  both  even  and  odd  extension  of  it’s  solution  u(x,  t )  into  ( -T ,  0)  we  see  that  the  operator 
changes  its  form  in  (-T,  0)  due  to  the  presence  of  the  first  derivative  with  respect  to  t.  Thus, 
parabolic  inverse  problems  are  usually  considered  in  {-T,  T )  instead  of  (0,  T)  [8],  [  13]-[  15],  [19]. 
Although  there  are  some  uniqueness  results  for  parabolic  inverse  problems  in  (0,  T),  but  they  are 
obtained  via  reducing  a  parabolic  inverse  problem  to  a  hyperbolic  one,  using  an  analog  of  the  inverse 
Laplace  transform,  see  [14],  [15]  and  subsection  3.3.1  in  [19]. 

Our  main  idea  is  that  it  is  possible  to  construct  a  proper  extension  for  the  non-stationary 
transport  equation  due  to  the  presence  of  the  parameter  v,  that  represents  the  vector  of  particle 
velocity.  To  illustrate  our  idea  consider  the  simplified  form  of  the  transport  equation 

u^{x,t,v)  +  (v,Vu+(x,t,v))  +  a(x,v)u+{x,t,v)  =  0,  (LI) 
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in  R"  x  (-7,0)  x  S'7,  S'7  =  {v  g  R"  :  Ivl  =  1}.  Consider  the  following  extension 
u~(x,t,v)  =  u+(x,-t,-v )  V(jc, t,v)  g  R"  x  (-7,0)  x  Sn. 

The  equation  (1.1)  becomes 

uj(x,t,v )  +  (y,Vu~(x,t,v))  -  a(x,  -v)u~(x,  t,v)  =  0,  in  R"  x  (-7, 0)  x  S'7. 

Hence,  the  principal  part  of  the  differential  operator  stays  the  same. 

When  trying  to  obtain  the  Lipschitz  stability  for  the  Inverse  Problem  formulated  below,  the 
authors  have  discovered  that  quite  cumbersome  conditions  need  to  be  imposed.  Thus,  they  have 
decided  to  limit  the  scope  of  this  paper  to  the  topic  of  uniqueness  only.  In  section  2  the  main  result 
(Theorem  1)  is  formulated,  and  it  is  proven  in  section  3. 

Statement  of  the  main  result 


Denote 


Q  =  {i  e  R'7  :  Ixl  <  R},  S"  =  {v  e  R"  :  Ivl  =  1}, 

H  =  Q  x  5"  x  (-7,  7),  r  =  5Q  x  S'7  x  (-7,  T),  Z  =  Q  x  S'7, 


H+  =  QxSn  x  (0,  T),  r+  =  dQ.  x  S'1  x  (0,  T), 
H ~  =  Ox  Sn  x  (-7,0),  r-  =  5Q  x  5"  x  (-7,0), 

so  that 


H  =  H+  U  Hr  and  T  =  T+  U  r  . 


Also,  denote 

C{H)  =  {s  g  Ck(H)  :  D"ts(x,t,v)  g  C(H),  \a\  <  k } 

The  transport  equation  in  the  domain  H+  has  the  form  [6] 

ut  +  (v, Vn)  +  a(x,v)u  +j"  g(x,t,v,/i)u(x,t,n)da ^  =  F(x,t,v),  (2.1) 

S" 

rs*2  _ 

where  v  g  S"  is  a  unit  vector  of  the  particle  velocity,  u(x,  t,  v )  g  C  ( H+ )  is  the  density  of  particle 
flow,  a(x,  v)  is  an  absorption  coefficient,  7(x,  t,  v)  is  the  angular  density  of  sources,  g(x,  t,  v,  pi)  is  a 
scattering  indicatrix.,  and  (v,  Vn)  denotes  the  scalar  product  of  two  vectors  v  and  Vm. 

Consider  the  following  boundary  condition 

u(x,t,v)  =  p(x,t,v),  (2.2) 


for  (x,  t,v)  g  {(x,  ?,  v)  :  Ixl  =  R,  t  g  (0,7),  v  g  5",  (n(x),  v)  <  0}. 
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Here  (n(x),v)  is  the  scalar  product  of  the  outer  unit  normal  vector  n{x)  to  the  surface  6Q  and  the 
direction  of  the  velocity  v.  Hence,  only  incoming  radiation  is  given  at  the  boundary  in  this  case. 
Equation  (2.1)  with  the  boundary  condition  (2.2)  and  the  initial  condition  at  t  =  0 

u(x,  0,v)  =  /(x,v),  V(x,v)  e  Z,  (2.3) 

form  the  classical  forward  problem  for  the  transport  equation.  Uniqueness,  existence  and  stability 
results  for  this  problem  are  well  known,  see,  e.  g.,  Prilepko  and  Ivankov  [22]. 

Suppose  now  that  the  absorption  coefficient  a(x,  v)  is  unknown,  but  the  following  additional 
boundary  condition  is  given: 


u(x,t,v )  =  q{x,t,v), 


for  (x,  t,  v)  e  {(x,  t,v)  :  Ixl  =  R,  t  e  (0,1),  v  e  S'\  (n(x),v)  >  0}. 

The  function  q{x,  t,  v )  describes  the  outgoing  radiation  at  the  boundary. 

Remark  1.  When  proving  the  uniqueness  Theorem  1  for  our  inverse  problem,  we  naturally 
assume  the  existence  of  two  solutions  ci\ (x,  v)  and  ai{x,  v)  of  this  problem  and  then  show  that 
a  i  (x,  v)  =  <22 (x,  v).  Therefore,  we  also  assume  the  existence  of  corresponding  solutions  u  \  (x,  t,  v) 
and  uo(x,t,v)  of  the  forward  problem,  which,  in  particular  satisfy  the  natural  compatibility 
conditions  between  boundary  and  initial  data. 


Introduce  the  function  y(x,  t,  v) 


7(x,  t,  v) 


p(x,t,v),  if  (n(x),v)  <  0, 

q{x,t,v),  if  (n(x),v)  >  0. 


Hence, 


nlr+  =  y(x,  t,v),  V(x,  t,v)  e  dD.  x  (0,7)  x  5".  (2.4) 

Thus,  we  obtain  the  following  coefficient  inverse  problem  for  the  non-stationary  transport  equation. 

Inverse  Problem.  Given  the  initial  condition  (2.3)  and  the  lateral  boundary  data  (2.4),  determine 
the  coefficient  a(x,v)  of  the  equation  (2.1). 

Theorem  1  is  the  main  result  of  this  paper. 

Theorem  1.  Suppose  that  derivative  dtg  exists  in  H+  x  5"  and  Wcf  g\\ C{jpyS!l)  <  r\  for  k  =  0, 1, 
where  r\  is  a  positive  constant.  Let  \f{x,v)\  >  >'2,  where  r2  =  const  >  0.  Assume  that  there  exist  two 
pairs  of  functions  (a\,u\ )  and  {ai,  M2)  satisfying  (2.1),  (2.3),  (2.4)  and  such  that 

<2i,<22  g  C(Z)  and  Ui,Uit,Uitt,Vui,S/uit  e  C(H+),i  =  1,2. 


Suppose  also  that 

[(af)(x,v)]2  =  [{af){x-v)f,i  =  1,2.  (2.5) 

Let  ||  uu  ||  c(TT)  -  r3,  where  r 3  =  const  >  0.  Then  there  exists  a  number  To  =  To(R,  tx,  )  >  R  such 
that  ifT  >  To,  then  a\  =  a  2  in  Zand  u\  =  U2  in  H+.  If  the  function  U2(x,  t,  1)  =£  0  in  H+,  then  it  is 

sufficient  to  have  T  >  R. 

Remark  2.  In  particular,  the  condition  (2.5)  is  satisfied  if  the  functions  a,(x,v),  (/  =  1,2)  and 
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f(x,v )  satisfy  one  of  the  following  pairs  of  relations: 


afx,v)  =  dj(x,  —  v)  and  /(x,v)  =  f(x,-v), 


or 


a i (x,  v )  -  a, (x,  -v)  and  f(x,  v)  =  -/(x,  -v ) , 


or 


a,(x,v)  =  -a,(x,-v)  and  /(x,v)  =/(x,-v), 


or 


a , (x,  v )  =  -a,(x, -v)  and  f(x,v)  =  -f(x,-v). 

Remark  3.  The  uniqueness  for  the  Inverse  Problem  considered  for  t  e  ( —T ,  T)  with  initial 
condition  given  at  t  =  0,  though  obtained  as  a  corollary  of  the  Lipschitz  stability  in  [18],  was 
obtained  under  stronger  requirements  for  the  differentiability  of  the  coefficients  of  the  transport 
equation  (2.1),  initial  and  boundary  data,  then  the  requirements  of  the  Theorem  1.  Note,  that  the 
proof  of  the  Theorem  1  is  valid  for  the  case  when  t  e  (-T,  T)  and  the  initial  data  is  given  at  t  =  0, 
thus  the  assumptions  of  the  Theorem  1,  being  naturally  extended  to  the  case  of  t  e  (-7',  T),  are 
sufficient  for  the  global  uniqueness  of  the  Inverse  Problem  considered  for  t  e  (-7',  T). 

Let 


Lou  =  u,  +  (v,Vn)  =  ut  +y~'J  ViUi, 
where  w,  =  du/dxi.  Introduce  the  function 


7=1 


t//(x,  t)  =  Ixl2  -  r]t2,  rj  =  const  e  (0, 1). 

Let  c  =  const  >  0.  Denote 

Gc  =  {(x,  t)  :  Ixl2  -  r)t 2  >  c2  and  Ixl  <  7?}. 
Introduce  the  Carleman  Weight  Function  (CWF)  as 

C  (x,t)  =  exp[Ai//(x,  t)]. 


Lemma  1.  Choose  the  number  ij  such  that  i)  e  (0, 1)  and  T  >  R/  Jrf .  Also,  choose  the  constant 
c  >  0  such  that  Gc  c  Q  x  (-7)  T).  Then  there  exist  positive  constants  Ao  =  Ao(Gr)  and 
M  =  M(GC),  depending  only  on  the  domain  Gc,  such  that  the  following  pointwise  Carleman 
estimate  holds  in  Gc  x  5"  for  all  functions  u(x,  t,  v)  e  C1  (Gc)  x  C(S")  and  for  all  A  >  Ao(Gc) 

(Lqu)2C2  >  2A(1  -  rj)u2C2  +  V  •  U  +  Vt, 
where  the  vector  function  (U,  V )  satisfies  the  estimate 
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\{U,V)\  <  MforC2, 


(2.6) 


and 


V(x,t,v)  =  2\(rjt  -  (v,x))w2C2.  (2.7) 

The  proof  of  this  lemma  can  be  found  in  [17].  The  following  lemma  provides  an  estimate  from 

the  above  for  an  integral  containing  the  CWF. 

Lemma  2.  For  all  functions  s  g  C(Gc)  and  for  all  a  >  1 ,  the  following  estimate  holds 


|  s(x,  t )dz 


g,  L  o 

See  Section  3.1  in  [19]  for  the  proof. 


C2 (x, t)dxdt  <  •  |  (s2C2)(x, t)dxdt. 


Lemma  3.  Let  T  >  R.  Then  for  any  c  g  (0  ,R)  there  exists  a  r/o  =  ijo(R,  T,c)  g  (0, 1)  such  that 


Gf  cQx  (~T,T)  for  all  rj  g  (rio(R,  T,c),  1). 

Proof.  By  the  definition  of  the  set  Gc 

Gcc  {fix  (~T,T)}  «max  y/(x,T)  <  c2-, 

8Q 

i.e.  when 


R2  -  r)T 2  <  c2. 


which  leads  to  the  following  inequality 

h  >  no 

Since  c  g  (0 ,R)  and  R  <  T  then  770  g  (0, 1).D 

Proof  of  Theorem  1 


R2-c2 

T2 


Suppose  that  there  exist  two  solutions  of  the  Inverse  Problem,  (ai,m)  and  (<22,  ui)-  Denote 

a  =  a\  -ao  and  u  =  u\  -112.  (3.1) 

From  relations  (2.1),  (2.3),  (2.4)  and  (3.1),  noticing  that  a\u\  -  a2U2  =  a{u  +  du2,  we  obtain 

ut  +  (v, Vn)  +  ai(x,v)u  +J"  g(x,t,v,/j.)u(x,t,/j,)doII  =  -au2,  (3.2) 

S" 

u{x, 0,v)  =  0,  V(x, v)  g  Z,  (3.3) 
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zdr+  —  0. 


(3.4) 


Denote 


u+(x,  t,v)  =  u(x,  t,v),  in  H+, 

Hence,  (3.2)-(3.4)  become 

uf  +  (v, Vz?+)  +  ci\(x, v)u+  +|  g(x,t,v,n)u+(x,t,n)dop  =  -du2, 


(3.5) 


u+(x,  0,v)  =  0,  V(x,  v)  e  Z, 


U  +  \r+  =  0. 


(3.6) 


To  apply  the  Carleman  estimate  of  Lemma  1,  we  need  to  construct  such  an  extension  u~(x,t,v) 
of  the  function  w+(x,  t,  v )  into  the  domain  H  ,  that  the  principal  part  of  the  differential  operator  of  the 
equation  for  the  function  u  (x,  t,  v )  in  H~  would  be  the  same  as  the  one  for  u  (x,  t,  v). 

Set 


Hence, 


Also,  denote 


and 


u~{x,t,v)  :=  u+(x,-t,-v)  V (x, t, v)  e  H  . 

uj(x,t,v)  =  -U((x,-t,-v)  V (x,  t,  v)  e  H  . 

ut(x,t,v)  :=  u2{x,t,v)  V (x, t, v)  e  //+ 


(3.7a) 


(3.7b) 


u2(x,t,v)  :=  u2(x,-t,-v)  V(x, t, v)  e  //  . 

In  the  integral  term  in  (3.2)  change  variables  /i  ->  -//  and  rewrite  (3.2)  in  the  form 

u'l(x,t,v)  +  (v, Vn+(x, t,v))  +  ai(x,v)n+(x,t,v)  +|  g(x,t,v,-^)u+(x,t,~ij.)dat 


=  -a(x,  v )  w  2  (x,  t,  v). 

Next,  replace  in  this  equation  v  with  -v.  We  obtain 

uf(x,t,-v)  -  (v,Vu+(x,t,-v))  +  ai(x,-v)u+(x,t,-v)  +j"  g(x, t, -v, -/./)«  '(x, t, -u)daL 


sn 


=  -a(x, -v)u2(x,  t, -v). 
Now  substitute  -t  for  t  assuming  that  t  <  0, 
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ut (x, -t, -v)  -  (v,  Vn+(.x,  -t, -v))  +  a i (x,  -v)n+  (x, -t, -v)  + J  g(x, -t, -v, -fu)u+(x, -t, ~/j,)doh 


=  -a{x,-v)uo{x,-t,-v). 


By  (3.7a, b)  we  obtain  in  H 


u t  (x,  t,  v)  +  (v,Vn  (x,  t,  v))  -  ai  (x, -v)n  (x,  t,  v)  - 


(3.8) 


-j"  g(x,-t,-v,~n)u  (x, t, n)cla u  =  a(x,-v)u2(x,-t,-v). 

S" 


Although  the  non-principal  part  of  the  operator  for  the  function  u~  is  different  from  the  one  for 
the  function  u+,  but  this  doesn't  affect  the  method  of  Carleman  estimates.  This  is  because  the 
Carleman  estimate  depends  only  on  the  principal  part  of  the  differential  operator.  Also,  by  (3.6), 


u  (x,0,v)  =  m+(.x,0,-v)  =  0  V(.x,v)  e  0x5". 

Thus,  we  have  constructed  a  continuous  function 


u(x,  t,  v) 


u+(x,t,v)  in  H+, 
u~(x,t,v)  in  H  . 


(3.9) 


Also,  note  that 


n(.x,0,v)  =  0  and  nip  =  0. 


Since  0  <  n  <  ]/(x,  v) I  in  Z,  there  exists  a  small  number  e  g  (0 ,R],  such  that 


uo(x,t,v)  ^  0,  V(.x, t,v)  g  H^, 


where  HI  =  Ox  (0,  e)  x  5". 
Thus,  we  can  rewrite  (3.5)  as 


a  =  u2 


it  +  (v, Vn+)  +  a\{x, v)u+  +j"  g{x,t,v, n)u+(x,t, fi)da 


sn 


(3.10) 


in  Ht.  Since  the  function  a(x,v)  is  independent  on  t,  we  can  eliminate  a  from  the  equation  (3.10)  by 
differentiating  this  equation  with  respect  to  t.  Thus,  denoting 


k+(x,t,v)  =  i^t-^j{x,t,v), 


(3.11) 


we  obtain 
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3.12 


+  (v,Vut)  +  ai(x,v)ut  +J  (gtu+  +gtut)dotl  = 


k+(x,t,  v) 


it  +  (v,Vm+)  +  ci\{x,v)u+  +|  g(x,t,v, gi)u+ {x,t, n)do h 


sn 


Denote 


w+(x,t,v)  =  ut  -  k+u+,  in  Ht . 


Hence  (3.12)  becomes 


wt  +  (v,  Vw+)  =  -ktu+  -  (v,Vk+)u+  -  al(x,v)w+  - 


Since 


where 


-j"  gw+do u  +  k+  |  gu+dot,  -j"  g,u+dGfl  +j"  gu+dot,. 

S"  sn  S"  sn 


u+(x,t,v)  N(x,t,T,v)w+(x,z,v)dz, 
0 


N+{x,t,  T,v) 


ut{x,t,v) 
ut(x,T,v)  ' 


Turning  equation  (3.14)  into  inequality,  we  obtain 

I  wt  +  (v,  Vw+)l  < 


(3.13) 

(3.14) 


(3.15) 


(3.16) 


<  K 


\w+\  +|  Iw+(.x,t,  v)\dz  +| 

0  S" 


\w+\do^  +||  \w+(x,T,v)\dTdot 

S"  o 


Let 


c2  =  R2  -  77s2, 


(3.17) 


consider  the  set  GJ ,  where 

Gt  :=  gc  n  {t  >  0}  =  {ixi2  -  rjt2  >  r2  -  Tje 2y  n  {(x,  t)  ■.  ixi  <  r,  t  >  o}, 


(see  Fig.  1.). 
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Fig.  1 .  1  -D  schematic  representation  of  the  domains  G+  and  G, . 


Also,  denote  P+  =  GJ  x  S".  Clearly, 

P+c  cHJ. 


Multiplying  both  parts  of  the  inequality  (3.16)  by  the  CWF,  squaring  both  sides  and  integrating  over 
Pc,  we  obtain 

|  \w,  +  (v,Vw+)\2C2dh  < 

Pi 


<  K 


Pi  i— 


— i  2 


I W+ 1  + 


|  l>v+(x, T,v)\dr  +|  \w+\do^  +j*j*  \w+(x,t, v)\drdat 


S"  o 


C  2dh, 


where  dh  =  dxdavdt.  Here  and  below  in  this  proof  K  denotes  different  positive  constants  depending 
on  numbers  r  \ ,  r2 ,  r 3 ,  s  and  norms 


al  II  C(Z)  5  llM2f  II  C(ff+)  ’  ll^M2  ||  c(ff+)- 


We  obtain 
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I  \wt +  (v,Vw+)\2C2dh  < 


(3.18) 


<K  \  \  \w+\ 


r 

+  I  \w+(x,T,v)\dz  C~dh  + 


+  K  |  |  \w+\2dotl  +|  [  |  \w+(x,T,v)\dz  do  ^  C 


Noticing  that  for  any  function  s(x,  t,  v)  e  C{H+ ) 


J"  |  J"  s2dojl  C2J/i  =  A  s2C2dh, 
pt  \s"  J  Pi 


where  A  is  the  area  of  the  unit  sphere  S'1,  we  remove  the  inner  integrals  over  S"  in  (3.18).  Hence, 
(3.18)  becomes 

|  \w~$  +  (v, Vw+)\2C2dh  <  K  J"  lw+l2  +  j"  \w+  (x,r ,v)\dr  C 2dh. 


Applying  the  Carleman  estimate  of  Lemma  1,  we  obtain 


2A(1  -  77)  |  \w+\2C2dh  +|  [V»Ui  +  Vu\dh 


K  |  \w+\2C2dh  +  K  |  |  lvv+(x,  t,  v)\dr  C 2dh. 


where  functions  U 1  and  V\  are  functions  U  and  V  from  the  Carleman  estimate  of  Lemma  1,  applied 
to  the  operator 


wf  +  (v,  Vw+). 


Using  the  Gauss’  formula,  we  obtain 


2A(1  —  77)  |  \w+\2C2dh  +  j*  ((t/i, V\),n{x, t))dS 
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< 


K  |  \w+\2C2dh  +  K  |  (  |  I w+(x,T,v)\dT  C 2dh, 


Pi  \  o 


where  ((C/i,  V\),n(x))  is  the  scalar  product  of  vector  function  {U \,V\)  and  the  unit  outer  normal 
n(x,  t)  to  the  boundary  dP^.  Hence,  using  Lemma  2,  we  obtain  for  2  >  1 


22(1-17)  |  lw+l2C2r^/7  +  J*  ((U\,Vi),n(x,t))dS 


< 


dPi 


<  K  |  \w+\2C2dh. 


(3.19) 


Choose  2o  such  that  7C/(22o(  1  -  77))  <  1/2.  Then  for  all  2  >  max{2o,  1}  we  have  from  (3.19) 

2(1  —  77)  |  \w+\2C2dh  +  ^  ((U\,Vi),n(x,t))dS  <  0.  (3.20) 

P;  8Pi 

Consider  the  boundary  8G +  of  the  set  G+ .  Denote 

dGl+  =  {(x,0  :  Ixl  =R}nG+, 

dG2+  =  {(x,  t)  :  Ixl2  -  rjt 2  =  c2}  fl  Gc, 


Hence, 


dG2+  =  {(x, t)  :  t  =  0}  n  G+. 

dGt  =  dGl+  U  dG2+  U  dG2+. 


(3.21) 


Also,  denote 

8P?  =  dG?  x  Sn,  i  =1,2,3. 

Thus,  we  obtain  from  (3.4)  and  (3.20) 

2(1-77)  |  \w+\2C2dh  +  J"  ((t/i,  Ci),/7(x,  t))J5  -  J  Vi(x,0,v)dxdav  <  0. 

Pt  sfL 

Denote 

H)  =  Qx(-s,0)xS",  Gc  :=  Gc  fi  {t  <  0},  Pc  =  GcxSn. 

For  c  satisfying  (3.17) 


(3.22) 


Pc  C  //-. 


Also,  let 
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(3.23) 


dG\  =  {(x,  t)  :  Ixl  =  R}  f|  Gc, 


dG2  =  {(x,t)  :  Ixl2  -  rjt2  =  c2}  fl  Gc, 


Hence, 


Let 


dGl  =  {(x,0  :  t  =  0}  n  Gc. 
dGc  =  oG  'c  U  dG2-  U  dG3-. 


dP{  =  dG\.  x  S'\  i  =  1,2,3. 


Denote 


Then 


where 


k  (x,t,v)  =  (“^|-)(x, t,v). 


u  (x, t,v)  =|  N  (x, t, t, v)w  (x,T,v)dr, 
0 


N  (x,  7,  r,  v) 


u2(x,t,v) 
Uo(x,T,v)  ’ 


Considering  the  equation  (3.8)  in  the  set  H ,  and  using  the  same  reasoning  as  above,  we  obtain 


22(1-77)  |  \w  \2C2dh  +  |  ((C2,  V2),n{x, t))dS  +  j"  V2(x,0,v)dxdov  <  0,  (3.24) 


dPi 


where 


Hence 


w  (x,  t,v)  =  u,  -  k  u  ,  in  Hc. 


(3.25) 


m  (x,  t,  v)  =|  N  (x,t,z,v)w  (x,T,v)dr,  (3.26) 

0 


and  functions  Ui  and  V3  are  the  functions  U  and  V  from  the  Carleman  estimate  of  Lemma  1,  applied 
to  the  operator 


w,  +  (v,Vw  ). 

Consider  the  set  Pc  =  P~  U  Pt-  By  (3.21)  and  (3.23) 
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dPf-  =  dPl+  =  {x  :  c  <  Ixl  <  R}  X  5"  :=  SP2. 

Hence,  summing  up  the  inequalities  (3.22)  and  (3.24),  we  obtain 

2(1-77)  |  lwl2C2J/7  +  |  ((f/i, Vi),n(x,t))dS  +  j"  (((/2, V2),n(x,t))dS  + 

Pc  8Pp  8 Pp 


+  |  [V2(x,0,v)  -  Vi(x,0,v)]dxdov  <  0. 

dPl 

Here  the  function  w  is  defined  as 


w(x,  t,  v) 


w+(x,t,v)  in  H+, 
w~{x,t,v)  in  H  . 


(3.27) 


We  first  need  to  cancel  out  the  possibly  non-zero  integral 

|  [V2(x,0,v)  -  Vi(x,0,v)\dxd(jv.  (3.28) 

8Pc 

Actually,  this  is  the  central  point  of  the  proof  compared  with  the  case  of  considering  ( -T ,  T)  instead 
of  (0,  T)  (see  Introduction).  By  (2.7) 

Vi(x,0,v)  =  -22(v,x)(w+(x,0,v))2C2, 

and 

V2(x,0,v)  =  -22(v,x)(w_(x,0,  v))2C2. 

Hence,  if 


[(w+(x,0,  v)]2  =  \yv{x,  0,  v)]2, 

then 

V2(x,0,v)- Vi(.x,0,v)  =  0  in  dPc- 


By  (3.9),  (3.13)  and  (3.25) 

w+(x,0,v)  =  ut(x, 0,v)  and  w~(x,0,v)  =  uj(x, 0,v), 
and  by  (2.3),  (3.5)  and  (3.8) 

w+(x,  0,  v)  =  -a(x,v)w2(x,0,v)  =  -a(x,v)/(x,v), 

and 


(3.29) 
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(3.30) 


w  (x,0,v)  =  a(x,-v)u2(x,0,-v )  =  a (x,  -v)f{x,  -v) 
for  all  (x,  v)  g  P\.  Thus  (3.29)  is  equivalent  to 

[a(x,  v)/(x,  v)]2  =  [ a (x, -v )/(x, -v ) ] 2 ,  V(x,v)  g  P2. 
Because  of  (3.30),  the  integral  (3.28)  equals  zero  and  the  inequality  (3.27)  becomes 

A(  1  —  77)  |  \w\2C2dh  +  j"  ((t/i, V\),n{x, t))dS  + 

+  |  (( U2,V2),n(x,t))dS  <  0. 

dP~- 

The  inequality  (3.31)  leads  to 

2A(1  -77)  |  \w\2C2dh  <  |  \(Ui,Vi)\dS  + 

Pc  3P?+ 

+  {  \{U2,V2)\dS. 

8Pl- 

Denote  dP2c  =  dP^  U  dP2c+.  Hence,  using  the  estimate  (2.6),  we  obtain  from  (3.32) 

2A(1  -  77)  |  \w\2C2dh  <  KX  |  \w\2C2dS. 

Pc  8P2c 

Let  S  g  (0,  s)  be  an  arbitrary  number.  Denote 

Gcs  =  {Ixl2  -  T]t 2  >  R2  -  rjS2}  n  {Ixl  <  R }, 
and  Pcs  =  Gcs  x  Sn  (see  Fig.  2). 


(3.31) 


(3.32) 


(3.33) 
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Fig.  2.  Schematic  representation  of  the  sets  Gc  and  Gcs  in  the  1-D  case. 
By  (3.17)  PcS  a  Pc.  Hence, 


Since 


|  lvvl2C2d/i  <j"  lvvl2C2J/7. 

PcS  Pc 


C 2{x,t)  >  exp[2X(R2  -  rjd2)]  in  Pc§, 

and 

C2(x,  t)  =  exp[2X(R2  -  rj£2)]  on  dP2c, 


we  obtain  from  (3.33)  and  (3.34) 

exp[22.(R2  -  r/S2)]  J  \w\2dk  <  Kexp\2X(R2  -  77s2)]  J  \w\2dS. 

Pcs  dP2c 

Dividing  this  inequality  by  exp\2X{R2  -  T]S2)]  and  letting  X  ->  00,  we  obtain 

|  lrvl2<i/7  <  0, 

PcS 


and  thus 


(3.34) 
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w(x,t,v)  =  0  in  Pc5. 


Since  S  g  (0,  s)  is  arbitrary  number,  then 

w(x,t,v )  =  0  in  Pc. 

Thus,  by  (3.15)  and  (3.26)  we  have 

u(x,t,v)  =  0  in  Pc. 

Since  rj  g  (0, 1)  is  an  arbitrary  number,  then  (3.10),  (3.35)  and  (3.36)  imply  that 

a(x,  v)  =  0  in  {Jr2  -  £2  <  Ixl  <  R)  x  5". 

Let  to  g  (-T,  T)  be  a  number  which  we  will  choose  later.  Denote 

GE(to)  :=  {(x, ?)  :  Ixl2  -  ^(t-  to)2  >  R2  -  77s2,  I.tI  <  R}. 
We  need  to  have  G£(to)  c=  {t  :  Itl  <  T}  x  Q.  Since 

to  -  s  <  t  <  to  +  £  in  GE(to), 

we  choose  to  g  (-T  +  £,  T-  s).  Also,  denote  Pc(to )  :=  Gf(to)  x  5"  (see  Fig.  3). 


Fig.  3.  Schematic  representation  of  the  sets  Gc  and  Gr(fo)  in  1-D  case. 
By  (3.37),  the  equations  (3.5)  and  (3.8)  in  the  domain  Pc(to )  become  respectively 


(3.35) 

(3.36) 

(3.37) 
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it  +  (v,Vm+)  +  al(x,v)u+  +|  g(x,t,v, n)u+ {x,t, jj)do #  =  0, 


for  t  >  0,  and 


u,  (x,  t,  v)  +  (v,  Vu  (x,  t,  v))  -  Cl\  (x, -v)u  (x,  t,  v)  - 


g{x,-t,-v,-jj)u  (x, t, n)do ^  =  0, 

S" 

for  t  <  0.  Since  these  are  the  homogeneous  equations,  we  do  not  need  to  divide  by  u 2  (x,  t,  v)  (see 
(3.1 1))  to  prove  that  the  conditions  (3.6)  and  (3.7)  imply  that  their  solutions  are  zeros.  In  other 
words,  using  the  above  method  but  without  the  introduction  of  the  functions  w+  and  w~,  we  establish 
that 

u+{x,t,v)  =  0  in  Pc(to)  n  {t  >  0} 

and 

u~(x,t,v)  =  0  in  Pc(to)  n  {t  <  0}. 

Therefore, 

u(x,t,v )  =  0  in  {(x,  t)  :  Jr2  -  s2  <  Ixl  <  R,  t  e  (-T  +  s,  T-e)}  x  Sn. 

This  means  that  now  we  can  consider  equations  (3.5)  and  (3.8)  in  the  domains 

{(x,t)  :  Ixl  <  Jr2  -  e2 ,  t  e  (0,T-s)}  x  Sn, 

and 

{(x,t)  :  Ixl  <  Jr2  -  s2  ,  t  e  (-T+8,0)}  x  Sn, 
respectively,  (see  Fig.  4),  and  the  boundary  T  of  the  domain  H  should  be  replaced  with 

{ixl  =  Jr2  - s2 ,  t  g  (o,r-s)>  xS". 
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Fig.  4.  Schematic  representation  of  the  new  domains  for  1-D  case. 


Repeating  the  above  proof  in  this  new  domain,  we  conclude  that 

a(x,v)  =  0  in  {  Jr 2  -  2s2  <  Ixl  <  R}  x  Sn. 

Without  the  loss  of  generality  we  can  assume  that  R 2  =  ms2,  where  m  is  an  integer.  Repeating  the 
above  process  m  times,  we  obtain 


d(x,v)  =  OinZ.  (3.38) 

However,  to  make  sure  that  we  indeed  can  repeat  this  process  m  times,  we  need  to  have  the 
observation  interval  (0,  T)  to  be  sufficiently  large.  We  now  establish  an  estimate  from  below  for  the 
number  T.  We  need  to  have  uo(x,  t,v)  ^  0  in  HI .  We  have 

t 

U2(x,t,v)  =/(x,V)  +  |  U2t(x,  Z,v)dT. 

0 

Hence,  \u2(x,  t,  v)l  >  r2  -  tr 3  in  H+.  In  order  to  have  r2  -  tr 3  >  0,  we  need  to  take  t  >  r^r?,-  Hence, 
we  choose 

0  <  8  <  min^^-,Rj. 
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(3.39) 


To  repeat  that  process  m  times,  we  should  have 


T 


d2 

>  me  =  > 


R2 


which  gives  us  an  estimate  for  the  number  T  from  the  below.  In  particular,  if 

‘ R • 

then  we  can  choose  e  =  R  and  the  estimate  (3.39)  yields  T  >  R. 

Finally,  to  prove  that  u+(x,  t,  i’)  =  0  in  H+,  we  observe  that  because  of  (3.38)  the  equation  (3.5) 
becomes 


ut  +  (v, Vn+)  +  ai(x,v)u+  +j"  g(x,t,v, fi)u+(x,t, fi)do =  0. 

sn 

Using  the  standard  method  of  energy  estimates  with  the  initial  condition  (3.6)  and  the  zero 
boundary  condition  u+  |  p+  =  0,  we  obtain 

u+(x,t,v)  =  0  in  H+. 
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